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We study spin relaxation and decoherence in nanotube quantum dots caused by electron-lattice 
and spin-orbit interaction and predict striking effects induced by magnetic fields B. For particular 
values of B, destructive interference occurs resulting in ultralong spin relaxation times T\ exceeding 
tens of seconds. For small phonon frequencies ui, we find a 1/-\/uj spin-phonon noise spectrum - a 
dissipation channel for spins in quantum dots - which can reduce Ti by many orders of magnitude. 
We show that nanotubes exhibit zero-field level splitting caused by spin-orbit interaction. This 
enables an all-electrical and phase-coherent control of spin. 

I. INTRODUCTION 

Although semiconductor spintronics is a field with already a substantial history as well as with commercial 
applications P spintronics with carbon-based materials is a young research area with excellent perspectives. Only 
very recently, a pronounced gate-controlled magnetoresistance response in carbon nanotubes connected to ferromag- 
netic leads has been reportedP Furthermore, spin injection and detection in single-wall carbon nanotubes has been 
demonstrated using a four-terminal geometry) 3 -' The interest to implement spintronic devices with carbon materials 
such as carbon nanotubes^ or graphene^ is mainly driven by the desire to improve material properties, for instance, 
for the spin relaxation behavior in these materials (as compared to more standard semiconductors like GaAs). This is 
so because carbon is a comparably light atom, thus, spin-orbit interaction is typically weakP Additionally, it consists 
predominantly of 12 C, which has zero nuclear spin, thus, spin decoherence and relaxation caused by the hyperfme 
interaction of the electron spin with the surrounding nuclear spins is weak. The advantageous material properties of 
carbon also trigger a large interest to create spin qubits^ in such materials. 

Here, we provide quantitative calculations of spin relaxation and spin decoherence times and show that they are 
dominated by a combination of spin-orbit and electron-phonon interaction. It turns out that such spin-orbit induced 
effects get strongly enhanced in small-radius nanotubes due to the curvature of the lattice, and result in energy 
splittings that even exceed those occurring in GaAs nanostructures. 

The interplay of such enhanced spin-orbit interaction with the one-dimensional nature of nanotubes results in a 
complex behavior with an extremely wide range of relaxation rates which can be varied over many orders of magnitude 
by an external magnetic field applied along the tube axis. We show that interference effects can result in ultralong spin 
relaxation times exceeding tens of seconds. By contrast, we uncover that for nanotube quantum dots, a spin-phonon 
dissipation channel exists with a sub-Ohmic spectral function (oc 1/y/uj, see below) which results in decreasing spin 
relaxation times for decreasing spin level splitting w. Compared to standard quantum dots (such as GaAs or In As 
semiconductors) this is a most surprising behavior, since usually the spin decay times increase for decreasing j 8 l 9 l 10 l 

Most remarkably, at zero magnetic field, the spin-orbit interaction induces a zero-field splitting in the energy 
spectrum. We show that this opens the door for an all-electrical control of spin in nanotube quantum dots, again 
based on the strong spin-orbit interaction. This feature is most interesting for spintronics applications where one 
aims at a spin manipulation without making us e of magnetic fields . Since quantum dots in semiconducting carbon 
nanotubes have been realized by several groups,^ 1 ' 12 ' 13 ' 14 ' 15 ' 16 ' 17 ' 18 -! we believe that our predictions are well within 
experimental reach. 

The paper is organized as follows. In Section [TTJ we introduce a theoretical model for a nanotube quantum dot 
and solve the spectral problem of the Hamiltonian of such a system. In Section |III[ we study spin-orbit coupling in 
nanotubes, consider different contributions to the spin-orbit coupling, and investigate zero-field-level splitting induced 
by spin-orbit coupling. In Section [rv] electron-phonon coupling in nanotubes is considered. Analitical expressions for 
the coupling of an elecron to three deformational accoustic phonon modes are obtained. In Section [Vj spin relaxation 
of an electron in a nanotube quantum dot is investigated and discussed. 

II. THEORETICAL MODEL 

We consider a single wall nanotube (NT) defined by the chiral vector C/j = niai + 712^2, where ai = ao(l,0) 
and &2 = ao(l/2, -\/3/2) are the primitive lattice vectors {clq — 0.246 nm) and n±, n2 € ZP The indices (ni,n2) 
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determine the radius of a NT R = |C^|/27r = ao\Jn\ + n\ + nin2/2n and the chiral angle (direction angle of 
C/t) # = arctan[V3ri2/(2ni + 712)] (see Fig. [TJ. Neglecting curvature effects (which lead to an inessential shift 
of the valley minima in k-space^) and SOI, we describe the system at the K = (27r/ao)(l/3, \/y/2) and the 
K' = (27r/flto)(-l/3, 1/a/3) point of the Brillouin zone (see Inset in Fig. [T]) by the Hamiltonian of grapheneP^l 

H = hv(T 3 k x a 1 + k y a 2 ) = hv(^^ K+ . k ^ elT3g ^ l ® e j, (1) 

where v is the Fermi velocity in a NT (v — 8.1 x 10 7 cm/sjP, (7j are Pauli matrices operating on sublattice space, 
and T3 = 1 (r 3 = —1) for the K (K 7 ) point, k is the electron wave-vector component along £ and, k is along (see 
Fig. [lj. It is convenient to perform a unitary transformation to remove the dependence on the chirality angle 9 from 
the Hamiltonian, i.e.: 

/ piraB ft \ 

u = ( e ;), (2) 

H = UHqU- 1 = hv(r 3 Ka 1 + ka 2 ). (3) 

Eigenvalues and eigenfunctions (in the rotated reference frame (<p,C)) °f the Hamiltonian ([3| at zero magnetic field 
are given by 

E K!k = ±hvy/ K * + k*, (4) 
p iK<,) - r / (') \ 



_ (K - tfe) , _ (K+Zfc) 

where r = (i?<p cos — £ sin 6>, i?<yS sin 9 + £ cos 0) . 

Periodic boundary conditions along the NT circumference [^(r + C^) = ^(r)] quantize the wave vector associated 
with the Ch direction [(k + K^) • = 27tto, m € Z]: k — > (m — r 3 i//3)/i?, where m <E Z and ^ = 0,±1 is 
determined by ni — «2 = 3iV + ^ (iV € Z). 4 A NT with ^ = [e.g. (n.,ra) armchair NT] has zero band gap and is 
called a metallic NT. Such a NT is not suitable to confine particles due to the Klein paradox in gapless structures.^ 
Therefore, semiconducting NTs {y = ±1) arc more favourable for quantum dot realizations, and we focus on this case 
in the following. An additional feature of semiconducting NTs with v = ±1 is that they allow us to avoid the problem 
of energy degeneracy at the K and K' points by applying an Aharonov - Bohm flux $ab = BttR 2 through the NT 
cross section?! Lifting the degeneracy is crucial for spin qubit realizations with controlled interqubit exchange.^ The 
Aharonov — Bohm flux leads to a shift of the quantum number to — > to + <5>ab/&o ($0 = hc/\e\ is the flux quantum) 
and to a Zceman splitting -E Km ,fc — * E Kmt k.s ( = E Km ^ + S^huiz, where luz — \e\gB /2moc, = ±1/2 is the spin 
projection on the NT axis. Therefore, the energy spectrum and wavefunction of an electron in a NT are given by 

E Km ,k,s<; = E Km! k + S c hw z , (7) 
*t,k,s c (<P>Q = *i'i>,C)K>, (8) 

where K m = (to + ^ab/'&o — T 3^/3)/i?, = Mff-B /2tooc, and \Sq) the spin part of the wave function. 

Now we consider a quantum dot (QD) which is made of a NT by the deposition of top gates on the N F 16 * 17 * ^ (see 
Fig. [2^,). The spacing between the gates L defines the length of a QD. We describe the confinement by the rectangular 
potential (see Fig. ^p) : 

V(C) -{ Vg ' C < or C > i, (g) 
V ^> ~ \ 0, < C < L. 1 1 

Recent experimental realizations of a NT QTj 19 l 20 12U provide clear evidence favouring the rectangular confinement in 
a QD, since Fabry - Perot interference observed in such experiments is a testimony for a NT QD with a well-defined 
length. Note that we consider the experimentally more accessible case, when the length of a NT QD L is much larger 
than its radius R (L 100 nm). For such QDs, the step-like potential drop happens on a length scale much larger 
than the lattice constant. Therefore, it does not introduce intervalley scattering. 

Straightforward calculations show that the bottom of the m-th subband of the NT spectrum under the top gates 
ftu|/« m | + V g divides the spectrum between the gates for this subband into two parts. Above the energy S«|/t m | + V g 
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FIG. 1: Two-dimensional hexagonal lattice. Here ai and a.2 are the primitive lattice vectors, Ch = niai + ri2a2 is the chiral 
vector (in this figure, we show the chiral vector with m = 5 and ni = 1), 9 is the chiral angle, £ is along the NT axis, and Rip 
is the azimuthal direction of a NT. In gray inset, the first Brillouin zone is depicted, where Y = (0, 0) is the center of the zone; 
K and K' are non equivalent points in the Brillouin zone. 



(dashed blue line in Fig. [2J;) the spectrum is continuous E Kmik (\k\ > (\V g \/hv)^/X + 2hv\n m \/V g ) and below there is 
a discrete spectrum E Kmt k n : 

* Km ,k, S< - I EKm k Sc: (| fc | >kc ,ke R), {W) 

where k c — (\V g \/hv)^l + 2hv\K m \/V g and allowed values of the quantized wave vector fc„ along the NT axis are 
found from the transcendental equation 

, , T (Hv) 2 k n k n 

t&nk n L= — — — . 11 

E Km , kn (E Km , kn - Vg) - (hvYK l m 

Here, k„ = y/t^, — {E Km ^ n — V g ) 2 /(hv) 2 . The wavefunction of an electron in a NT QD can be written as follows 
where 

*m,fc(0, C < 0, 

*m,*(C) = <( *£,k(0, < C < (13) 



Here, for a discrete spectrum {k = fe„ < fc c 



< fcn (C) = Ae^f V (14) 
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FIG. 2: Nanotube quantum dot. a. Nanotube with two top gates, b. Longitudinal confinement potential, c. Scheme of the 
band structure. 



< fcn (C) = 



1 



(15) 
(16) 



where 



z 



(') 



-, _ . -mi k n K m ,-ik n n _ .« *Km, — »fc» ^l^miK 

z (>) _ z (0 ' _ z (') 

K 771 1 ^71 ^ 777, 7 ^ 77, ^ 771 1 ^ 7~l ^ 771 7 77, 



(') 



z 



(') 



2; 



(') 



B = Ale 



/■£ Tin . — trCj-y ^rn 



v(') 



(17) 
(18) 



and A can be found from the normalization condition 



1 = || * || a = M a [(*W _, s ) 2 + l]-^ + |B| 2 [(zW + 1 +4L|C| S 



+Rc 



Ktti , — ik n ' 

Sj) 



CD*[{zZ,- kn Y + A W {l-e 

ttv n 



-2ik n L 



For the K-point, we obtain 



C = a]- + ilxaC, D = A 1 -- tbnC, 



ImC = 



2 k n 



-i 



,k n K"m ™ri 



AE, 



E, 



K>m fori 



(19) 



(20) 
(21) 
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B = Acos(k n L) 



E Km , kn { E K m .k n - v g ) - (hvK m y 



For a continuous spectrum (|fc| > k c ), we make the following ansatz: 

<*(C) = ^ ( \ k ) +Re-*t ( Z ^f k 



< fc (C) = A c e^ [ Z Y 



where |fc| > (|y 9 |/fiu) V 1 + 2fiu|/s m |/V fl , fc = ±^[(£ Km , fc - %) /H 2 - Kn) 2 , and find that 

{z Km ,- k ~ z Km - k )(z Km - k - z Km _~ k )^K 2 rn + fcV K m + fc2 



kkcos(kL) — ism(kL)(y + k 2 ^n 2 n + k 2 — n 2 t ) 

Z K m ,k — Z Km ,-k 



,-k 



z K m ,-k z Krn ,-k z K m ,-k z Km -k 

R = A c + B c -1, T = A c e lkL + B c e- lkL . 



(22) 

(23) 
(24) 
(25) 



(26) 

(27) 
(28) 



III. SPIN-ORBIT INTERACTION IN NANOTUBES 

Next, we take spin-orbit interaction (SOI) effects into account. In graphene, there are two main mechanisms of SOI: 
Intrinsic SOI, H™q = Ai nt T3<73S z }^ and extrinsic SOI (Bychkov - Rashba like), H$q = (Ae + A curv )(T3<TiS y — (J2S X ) 
which is due to the asymmetric confinement potential normal to the graphene sheet (A^pSI anc j curvature induced 
effective electric field of rippled graphene (A curv )l23 (sj is the Pauli spin matrix). In a NT, i.e. a graphene sheet rolled 
up into a cylinder, the spin components perpendicular to the NT axis become dependent on the polar angle 



i (-S+e^ + S-e~ iv ) 
2S C , 



= S 4 



s. 



(29) 
(30) 
(31) 



where, in the eigenbasis of S ?J 2S ( \ f} - | T), 2S C | j) = -| |>, S+\ T) = S-\ I) = 0, S+\ j) - | T), and 5_| T> = | j). 
Therefore, for a NT, the intrinsic SOI Hamiltonian is given by 



flgg = A int T 3 a 3 (5 + e^ + S-e-^), 
the extrinsic SOI term due to Ae is given by 

Hi = A E [2T 3 a 1 S < -ia 2 (~S + e lv + S-e- 
and the extrinsic SOI term due to curvature of a NT is given bjBH 

= iA^ niv a 2 (S+e* + S-e~^) + Ajj^rs^, 



')] 



(32) 
(33) 
(34) 



where A^ = -A(V° p - V£,)a /8V3B£ vir and ALv = A(3t£ + 5V^ p )a /8V^Re^ (A = UmeV^e^ = 7.3 eV, 
V° p = 6.38 eV, and V^ p = -2.66 eV [3D]). Note that at moderate electric fields (E < 0.1 V/nm), the last SOI term is 
dominant (A int w 1 //c<Y.' J7 A E < A int , and A^ urv : 
be safely neglected. 

The last term cx aLivS'c ( wnere ALrv ~ 0.17meV/i?[nm]) in Eq. (34 1 leads to a shift k 



(0.26 meV/i?[nm])) and, therefore, the other types of SOI can 



Km i A c 



(where ± corresponds here to | f) and | J.) states) and, therefore, to a spin splitting: 



^xr+ ,fc„, + l/2 ^V-,fc„,-l/2 



Swz — 2sgn(m — r 3 ^/3)A|! 



(35) 
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(for \ > *n). Thus, SOI ex Ac Ur v acts as an effective magnetic field resulting in a level splitting (2Ac Urv ) at zero 
magnetic field, as has been now experimentally confirmed nowP^ Note that this zero-field splitting does not violate 
Kramers theorem, since time reversed states correspond to different non-equivalent K-points and are degenerate at 
zero -B-fields (see Fig.[3|. The existence of the zero-field splitting opens up an intriguing possibility for spin resonance 
experiments without any magnetic fields: the first term in Eq. ( |34[ ) allows electric-dipole transitions between spin- 
up and spin-down states, the second term (as an effective magnetic field) splits these states, and thus oscillating 
electric fields perpendicular to a NT lead to electric-dipole spin resonance with resonance frequency to = 2A| urv /7i ~ 
33 x 10 10 s" 1 and Rabi frequency to n w 1.6 x 10 5 s" 1 at E = 10 V/ cm and V g = 2.3 meV (see APPENDIX [A} . 




-3-2-1 1 2 3 

B (T) 



FIG. 3: Lowest energy levels of electrons in a NT QD at low magnetic fields (R w 1.6 nm, L = 100 nm, g = 2, and 
V g = fiv/iOR « 8.5 meV). Solid curves correspond to t he case of zero intervalley mixing. At zero magnetic fields there is 



splitting of the levels due to the second term in Eq. (341. The magnitude of the splitting is 2Ac Urv = 0.22 meV. The arrows 
indicate crossings of j f) and j J.) states of a certain K point. At these levelcrossings, a singularity appears in 1/Ti (see below 
in Fig. |6j. Dashed curves correspond to the case of weak intervalley mixing (A K -k' = 0.05 meV). It modifies the zero-field 

splitting 2|Acurv| — > 2y A|j urv + A|,_ K , and opens up avoided crossings (with the value 2|A K -k'|) of the levels with the 
same spin orientation but different valley index. This is illustrated in the inset which is a blow-up of the center region of the 
spectrum. 

We know that intervalley mixing leads to splitting of the levels corresponding to different K points which has 
been observed in Ref. |331 Such mixing does not split Kramers doublets (in the case of time reversal symmetric 



intervalley scattering) but modifies the magnitude of the splitting (2Ac Urv — > 2^/Ac Urv + A^ K ,, where Ak-K' is 
the intervalley mixing strength) between spin-up and spin-down states of a certain K-point and leads to anticrossings 
at non-zero B-ficlds of the levels with the same spin orientation but belonging to different K valleys (see Fig. [3] and 
APPENDIX p. 

In the case of negative voltage applied to the top gates (V g < 0), hole states become localized instead of electrons. It 
can be shown that the energy spectrum of the lowest levels of holes has the same structure as for electrons (illustrated 
in Fig. [3]) but shifted down by the energy gap E g w 2hvy/ + k2 w 228 meV. From Fig.[5]we see that electron energy 



whereas there are no crossings of the two highest 
ue to zero-field splitting of energy levels) breaks the 



levels cross at hioz = 2t3V>Acutv (indicated by arrows in Fig 
levels of holes at nonzero magnetic fields. Therefore, SOI (c 
electron-hole symmetry. For the estimation of the SOI constants we use band parameters of bulk graphite. Note that 
for small radius NT, due to curvature effects, strong hybridization of bands can modify the band paramenters of a NT 
and, thus, the SOI constants or the g-factor. If the SOI constant had the opposite sign due to hybridization, then the 
energy spectrum for electrons would look like the one for holes and vice versaP^. Hence, in the case of negative Acurv, 
there are crossings of levels for holes and not for electrons (at B ^ 0) . Such electron-hole asymmetry in the spectrum 
can provide us with information about the sign and the magnitude of the SOI constant and about the (/-factor. 
Now, we turn to the quantitative discussion of the spin relaxation time in nanotube quantum dots. We take the 
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first term in Eq. ( 34 1 into account in the framework of perturbation theory, which leads to the solution of the Dirac 
(eigenvalue) equation for the lowest levels (Hq ± huiz/2 + #go rv )V'Q,o,±i/2 = £o,o,±i/2'0Q,o,±i/2 m first order in H§q v : 

£o,0,±l/2 ~ ^±,fco,±l/2i (36) 

L f ±ao 

^b,o,±i/a(v.C) ~ ^±,^,±1/2(^0 + ^^*^,^1/2(^0 + ^ J k dk\^^ i k Tl/2 {^X), (37) 



, fc(OM<^ f , fco (0) , 

A± = iiA^— ^ ° . (38) 

Note that the function fe , (0) cr 2 < i ) j C T & (C) i s either symmetric or antisymmetric with respect to inversion at 

( = L/2. Hence, there is a selection rule for SOI between quantized levels, namely, 

(*^A/AHaU,*»(0) « 1 - s § n K - l/3)sgn(m - l/3)(-l)"'+". (39) 

Thus, A+ ( = (A~ , = 0) for odd (even) n'. 

IV. ELECTRON-PHONON COUPLING IN NANOTUBES 

For definiteness, we consider only such (77,1,77.2) NTs that ni —772 = 3A^+1 (y = 1). Then the two states (V>o,o,±i/2) 

with the lowest energy (at ^ab/'&o > Ac Urv -R//w) of a NT QD belong to the K point with t 3 — I (see Fig. I3I. 
Phonon induced transitions (which become allowed due to SOI) between these states give the dominant contribution 
to spin relaxation of a single particle in a NT QD. Despite quite complicated phonon dispersion relations in NTs it 
is possible to find analytical expressions for the electron-phonon coupling in NTs in the case of low-energy phonons. 

The splitting between the £0,0, ±1/2 states is less than 1 meV at B < 10 T. Phonons with much higher energies are 
not favourable for transitions between these levels. The energy of the radial breathing mode is TuJrbm > 8.5 meV 
|31j for NTs with R < 1.5 nm, which excludes that mode (and all higher modes) from our analysis. Thus, only three 
acoustic phonon modes are important for spin-flip transitions between the lowest two levels: the twisting mode (TM), 
the stretching mode (SM) , and the bending mode (BM) P» To describe these modes we use a continuum modeP in 
which the equation of motion for the displacement u(r, t) — {um, uq, u r ) is given by 



where the force-constant tensor 



A = 



u(r,*)=Au(r,*), (40) 



V 2 „, cfV^+lrVL ^#V r I (41) 



R 2 <pC I ^ R 2 <P<f 

\ £LV7 C > ~ 2c * V7 _L_ 

\ fl V r R V C R 2 

is invariant under the group symmetry operations of a NT ) 4 * 32 * Here, q and Ct are the longitudinal and transverse 



phonon velocities, respectively (q = 20.9 km/s and c t = 12.3 km/s [36 ). Substituting the solution of Eq. (40 1 in the 
form u(r, t) — A a exp[i(rmp + qC — tot)] (q and to are the phonon wave vector and frequency, respectively, and a is the 
phonon mode) and keeping only leading terms in qR (qR <C 1), we get for TM phonons (m = 0): 



for SM phonons (m = 0): 
and for BM phonons (m = 1): 



Lo T = c t q, A T = A T (1,0,0) J (42) 
lu s = c s q, A s = A s (0, 1, -iqRrj) , (43) 
lu b = c s Rq 2 /V2, (44) 

B 

2 



A B ( . ir}(qR) 2 . n , ii{qR) 2 , 
A b = -£ » + , 1 - J ^r L - » (45) 
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where c s = 2(c t /ci)^/c 2 - cf, r] = (cf - 2c 2 t )/cf; Aj = v / h/2Mu; : j (M is the NT mass). We see that TM and SM show 
linear dispersion, whereas BM exhibits quadratic dispersion. Note that these results are only valid for long-wavelength 
phonons (qR <C 1 and uj < wrbm)- 

The clectron-phonon coupling is expressed by the operator 



el— ph 



Vi v 2 
v: v x 



where for the K-point 



Vi = gi(u vv 
_ 1 du v 



), v 2 = 

u C( = 



92 e 
du $ 



+ H.c, 

214,^ — 



He 
du v 



+ 2m vC ), 
1 du<; 
Hlhp' 



(46) 

(47) 
(48) 



gi ~ 30 eV is the deformation potential constant (which appears in diagonal elements of V^ph), and the off-diagonal 
coupling constant g 2 ~ 1.5 eV (which is caused by change in the bond- length between neighboring carbon atoms) 
Using Eqs. (|42|-((46l, we get for the TM: 



V? = 0, V 2 T = -g 2 A T qe 



3iSgi(«C-WT*) 



for the SM: 



and for the BM: 



2 2 

Vf = 2ig 1 A s q%e l(qC -~ ulst) , V 2 S = ig 2 A s q^e 3ie e 1 ^'^ , 
c i c t 



(49) 



(50) 



Vy = V2g 1 A B q 2 R^e l( ^-^ t \ V 2 B = g 2 A B q 2 R^e 3ie e l ^ + < 1 < : -^ t \ 



(51) 



Note that the electron-phonon coupling in nanotubes is very strong (for example, compare g± ~ 30 eV with a 
deformational acoustic coupling constant in GaAs So ~ 6.5 eV). Furthermore, the electron wave function is highly 
localized in the dot region (it decays exponentially outside the dot). Thus, the phonons in the contacts or substrate 
can be safely ignored for our purposes. Moreover, we neglect the effect of the substrate on the phonon modes. This 
is justified due to the relatively weak coupling between the substrate and the NT, very high stiffness and rigidity of a 
NT, and, last but not least, very small atomic displacement amplitudes in an acoustic phonon wave (which is a few 
percents of Angstroms only) . 



V. SPIN RELAXATION IN NANOTUBES 

We are now able to analyze spin-flip transitions between the lowest energy levels induced by long-wavelength 



phonons. Using Eq. (37 1, the matrix element of such a transition is given by 

M u = (^0,0,-l/2|Kl-ph|V , 0, 0,1/2) 

= E {^(^ -, fc o,-i/ 2 l^i- P h|^ =iA _ 1/2 ) + ka/2 \V e :- P ^^ MA/2 )} ■ (52) 

k 

Here the sum stands for summation over the discrete k n and integration over the continuous |fc| > 
(\V g \/Hv)^l + 2hv\K m \/V s 



a- 



FromEq. (l2j, (*« roi , fen , ± i /2 |e <m ^e^|* Km3)fcn ,, ±1 / 2 ) = ^> mu k n \e i "^ m3iKl )8 mi , m2+m3 . Therefore, only phonon 
modes with m 2 = 1 give non-zero contribution to spin-flip transitions (this is an additional reason why we do not 
need to consider higher phonon modes with m 2 > 1). Thus, only BM-phonons are responsible for the spin relaxation, 
whereas TM- and SM-phonons (with m 2 = 0) cannot flip the spin. 

In the framework of Bloch - Redfield theorypZI the spin relaxation time induced by BM-phonons is given by 



t l L dq{2N - + mu2 ¥i^m [s{q - qo) + s(q + qo)] 

2 5 / 4 7TL |2 



H 2 y/csRuQ 



{2N U]0 + l)\M UJ0 \\ (53) 
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e. 



,fe ,-i/2l/ S ~ l^z - 2r 3 Al! urv /S|, g = 



is 



2uj /c s R, and 7V W = [exp(/kj//c B T) - l]" 1 

0(A*) for SM 



where w = l^+, fco ,+i/ 2 

the Bose distribution function. Note that pure dephasing 1/TL = for BM phonons and l/T v 
and TM phonons, therefore, I/T2 = 1/2T\ + l/T v — l/2T\ in first-order perturbation theory. 

We used the Markov and the secular approximations in the derivation of Eq. (53). We can estimate the correlation 
time in the phonon bath to be r c ~ lps. Therefore, the Markov approximation (Ti 3> t c ) and the secular approxima- 
tion (luqTx 3> 1) are valid except for the energy regime close to the level crossing at ui = 0. Moreover, our estimations 
of the electron-phonon coupling are valid for phonons with the wavelength shorter than the full length of the NT 
Int- Therefore, in the case of a small splitting between spin-up and spin-down states (long wavelength phonons), the 
results are trustworthy for sufficiently long NTs (IntQo ^S> 1), for example, if the spin splitting is 1/ieV, then the NT 
length should be greater than 700 ran. 

We now study spin relaxation induced by low-frequency phonons (u>o ~ \u>z — 2r3A|j urv | — > 0). As shown above, 
such spin relaxation occurs near the level crossing indicated by arrows in Fig. |J One can show that |M Wo | 2 cx u>q 
and N Uo cx T/uj (at k^T 3> Huiq) for uj — > 0. Moreover, the density of states for one-dimensional phonon modes 
with quadratic dispersion, i.e. the bending modes responsible for spin relaxation, has a van Hove singularity at zero 
frequency. It goes like where ujq is the phonon frequency of the bending mode. This translates into the 

existence of a singularity in the noise spectral function J(luq) cx l/y/Uo which describes particle spin relaxation due 
to coupling to NT lattice vibrations via SOI and electron-phonon interaction. Therefore, 

1/Ti cx l/^Jo (54) 

at low ujq. To the best of our knowledge, this is the first system that exhibits a 1/y/uJo spin-phonon noise spectrum 
at low frequencies. Such a result (fast relaxation times at small splitting between spin-up and spin-down levels) is 
counter-intuitive in the light of the commonly expected long Ti time for NTs (due to the expected weak SOI) and 
compared to the usual behaviour of the spin relaxation time (1/Ti oc ll>% at low magnetic fields) for conventional GaAs 
QDsP 
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FIG. 4: Chirality dependence of the electron spin relaxation for NT QDs. The spin relaxation rate as a function of a field 
for a (40,0) zigzag NT (solid curve, 9 = 0), (-24,47) NT (dotted curve, 9 w tt/2), and (13,32) NT (dashed curve, 9 » tt/4). 
R w 1.6 nm, L = 100 nm, g = 2, V g = hv/AQR w 8.5 meV, T = 0.1 K. 



To better understand Eq. (54 1, we consider the spectral density of the electron-phonon correlation function 

/>oo 

J mk {uo) = / di(m|Ki- P h(0) |fc) (k| Voi-ph(t) \m)e wt , 



where the overbar denotes the ensemble average. We first analyze this expression for GaAs QDs and later on for NT 
QDs. 
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FIG. 5: Spin relaxation of a (40,0) zigzag NT. The dependence of the electron and hole spin relaxation rates 1/Ti on a 
parallel magnetic field B (lower horizontal axis for electrons and upper one for holes) due to SO coupling of the two lowest 
states (^ x ± kQ ±1 / 2 ) to higher states of the discrete spectrum ( , 5 j< =f k ±1/2) * s snown by the dashed curve and that due to SO 

coupling to states of the continuous spectrum =1=1/2) by the dotted curve. The total spin relaxation rate (solid curves) 

exhibits destructive interference at each odd crossing of the curves related to the abovementioned contributions (R w 1.6 nm, 
L = 100 nm, g = 2, T = 0.1 K, V g — ±hv/AQR ~ ±8.5 meV, where the upper/lower sign is for electrons/holes). The electron 
(blue solid curve) and the hole (red solid curve) spin relaxation time map onto each other by a shift along the magnetic field 
axis by AB — 2Al urv /gfiB- 



For the phonon-induced relaxation rate between levels split by the Zeeman term, we find 1/Xi oc JyzioJz) oc 
EqEiC-Nu; + 1 / 2 )l( 1 l^e iqr |0(/|Hso V l 2 )l 2(S ( w - w z), where A u is the electron-phonon coupling strength and u is the 
phonon frequency. Therefore, the corellation function defines the phonon-induced electron spin relaxation times. Let 
d be the single phonon degree of freedom (related to the dimensionality of the underlying lattice structure). Then, 
for GaAs semiconductor structures with linear in momentum H<?q v , we get ^ q — > J dqq^ 1 J df2 q , (l|e lqr |i) oc q (in 

dipolc approximation), (l\HgQ V \2) oc Uz> N u oc 1/u) (at k-gT ^> ftw). Taking into account that oc for the 

coupling between an electron and a piezoelectric phonon, we obtain Jia(wz) w z +1 m the case °f linear dispersion 
of a phonon (lu oc q). For deformational acoustic phonons, oc i/cJ, therefore, Jn[oJz) oc w^" 1 " 3 . Therefore, at low 
frequency, the spectral density function of the electron-phonon coupling is super-Ohmic ( J {to) oc lu" n > 1) even for 
all phonons in all dimensions. 

This is fundamentally different for the NT QDs discussed here: Since i?gQ V in a NT couples spin to the azimuthal 
degree of freedom (see Eq. ( 34 1 ) and the azimuthal component of the phonon wave vector is quantized (see Eq. (51)), 
we get (l|Aje lqr |Z) (1\H^q v |2) oc 1 + O(q). Thus, for deformation-acoustic phonons (A u oc y/uj) with quadratic 
dispersion (J2q ~ * / duj/^/uj), we obtain Ju(ujo) oc 1/y/uJo and recover Eq. (54 1. The noise spectral function J 12(1*}) 
describes particle spin dissipation due to coupling to NT lattice vibrations (via~SOI and electron-phonon interaction). 

As shown in Fig. [6j the magnetic-field dependence of the spin relaxation rate of a NT QD is exceptional in 
comparison to that of a conventional semiconducting QD. First, there is a singularity of the electron spin relaxation 

rate at luq — > (or at u>z — > 2r3Ac U rv) in contrast to the usual super-Ohmic behavior of 1/Ti in GaAs or InAs 
QDs (compare Fig. ^ with Fig. 1 :n !!). Remarkably, the position of this symmetric singularity gives us a direct 

II 

measurement of the SOI constant A" urv and valley index T3 of an electron in a NT. The singularity is at positive 
magnetic fields for the K point (T3 = 1) and at negative magnetic fields for the K' point (7-3 = —1). In Fig. IrJl where 

I — 

T3 = 1 and 2Acurv ~ 0.22 meV, the singularity is at B m 1.9 T. If the SOI constants and g factors are the same for 

both electrons and holes, then the electron and hole spin relaxation curves map onto each other by a shift along the 
magnetic field axis by AB — 2AjJ urv /<?/iB (compare the blue and the red curves in Fig. [6j . 

We have also studied the chirality dependence of the spin relaxation rate as a function of the magnetic field. 
Different chirality nanotubes show qualitatively similar spin relaxation properties. In Fig.|3J the spin relaxation time 
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FIG. 6: Spin relaxation of a (40,0) zigzag NT. The dependence of the electron and hole spin relaxation rates 1/Ti on a 
parallel magnetic field B (lower horizontal axis for electrons and upper one for holes) due to SOI of the two lowest states 
{^f x ± kg ±i/ 2 ) to higher states of the discrete spectrum k ^1/2) ^ s shown by the dashed curve and that due to SOI to 

states of the continuous spectrum x t t ^1/2) by the dotted curve. The total spin relaxation rate (solid curves) exhibits 
destructive interference at each odd crossing of the curves related to the abovementioned contributions (R « 1.6nm, L = lOOnm, 
g = 2, T = 0.1 K, V g = ±hv/l50R w ±2.3 meV (a), and V g = ±hv/3R « ±113 meV (b), where the upper/lower sign is for 
electrons/holes). The electron (blue solid curve) and the hole (red solid curve) spin relaxation time coincide by a shift along 
the magnetic field axis by AB — 2Ac Urv /<?A t B. In Fig. [6k (small V g ), the continuous spectrrum substantially influences 1/Ti. 
In contrast, in Fig. [6]d (large V g with many bound states in the dot), the discrete spectrum mainly determines the magnetic 
field dependence of 1/Ti. 



for NT QDs with different chirality but approximately the same NT radius is shown. From this figure we conclude 
that T\ depends on the chirality of a NT, although it has the same qualitative behavior as a function of a magnetic 
field. 

VI. INTERFERENCE EFFECTS IN SPIN RELAXATION 

We note here that the spin relaxation rate for flat GaAs QDs in in-plane magnetic fields is a monotonic function of 
B (up to about 14T),™ whereas, as shown in Fig. [6j it oscillates with B for NT QDs. The oscillations are caused by 
interference effects of two types: (i) interference of a phonon wave in a NT electron cavity bounded by the confining 
potential V(£) due to top gates (see Fig. [2]).; (ii) interference between various contributions to the spin-flip transitions. 
For clarity, we will now study these two types of interference phenomena separately. 



A. Interference of phonon waves 

To illustrate the first effect, we only consider one contibution to spin-flip transitions, namely, that due to the first 



term in Eq. ( 52 1 



Wo kg 



,fco,-l/2 



Ki-phl*, 



-1/2/ 



(55) 



Note that At = due to selection rules (see Eq. (39 1). The corresponding spin relaxation rate 1/T[ due to this term 



only is shown in Fig. [7] 

From Fig. [7] we see that l/T[ exibits oscillations as a function of the ratio between the NT QD length L and the 
phonon wavelength A p h: qL = 2n(L/\ p i 1 ). We attribute such oscillations to interferences of the phonon wave in a 
NT electron cavity bounded by the confining potential V(Q due to the top gates (see Eq. Q and Fig. [2]). Such an 
interference effect is reminiscent of a Fabry - Perot - type interference of a phonon wave where the electron levels 
in the dot play the role of a cavity. The coupling between the phonon wave and the cavity is described by the 
electron-phonon interaction V^i-ph- 
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FIG. 7: Interference phenomena in the spin relaxation rate due to the first contribution to the spin-flip transition (oc M' UQ ). 
Here, q is the phonon wave vector and L is the length of the NT QD (R ~ 1.6 nm and T = 0.1 K). 



At the minima in Fig. [7] the coupling between the electron cavity and the phonon waves becomes small. For an 
ideal cavity (with no loss), the matrix element of the spin-flip transition goes to zero at the minima. For instance, in 
the case of a rectangular hard wall potential, the squared modulus of the phonon-induced spin-flip transition is given 

by 



1 [2 TTX 






\VZ sm T 






( 8tt 2 sin(gL/2) 
\qL[^ - (qL)* 



(56) 



which is zero at qL = 4ir, 6tt, 8ir, .... Therefore, electron-phonon coupling is switched-off at these interference 
minima. In the case of a NT QD with a rectangular confining potential with finite barriers, however, due to the 
penetration of the electron wave function into classically forbidden region, the electron-phonon coupling is small but 
nonzero at the minima of the matrix element of the phonon-induced transition and the minima are shifted from those 
for an ideal cavity. As can be seen from Fig. [7] this shift and the minimal values of the electron-phonon coupling are 
more pronounced with increasing the barriers hight V g . 

Interference effects in a NT QD occur only for confinement with well-defined length (for all bound states) and are 
absent for soft potentials such as parabolic confinement. Note however that the rectangular potential seems to be 
a good approximation for the confinement in a gated NT QD, since Fabry - Perot interferences (for electrons) have 
been observed in such a system!^ 



B. Coherence of different contributions to spin-flip process 



In this subsection, we study interference effects due to various contributions to the spin-flip transitions described 
by Eq. (52 1. Let us consider the case of weak confinement with small V g — 2.3 meV. In this case, Eq. (52 1 can be 
rewritten as follows: 



= M+ + M- +M+ + M-, 
K = (At 



(57) 



(^+,fc 0> l/2l^-Phl^ +,fco,l/2> 

MI = \7(f v ^_ 1/2 |Ki-ph|«V iifcl ,_ 1/2 ) 



r+,fc,l/2 



r+,fc ,l/2/' 



13 



AC = 



■o ,«0> 



-l/2l^el-ph|*^ =ilfc _i/ 2 > 



(fe± = (|V fl |/ft«)^/l + 2fiv|*£il/^)- 

the two lowest levels (£0.0, ±1/2) an( i higher discrete levels {£ 



Here and AC 



are contributions to the spin-flip transitions due to SOI of 



1,0. +1/2 an d £—1,1 -1/2)- Note that the coupling to other 
higher levels is forbidden by the selection rule Eq. ( 39 1 . The contribution of these two terms to the spin relaxation 
rate is shown in Fig. [8j It can be seen that these two terms interfere (constructively) which leads to a change in 
the amplitude and period of the oscillations. Such constructive interference (see Figj8J) between M~£ and AC just 
increases the spin relaxation rate and, therefore, is not the dominant one. But next we consider a different interference 
effect which reduces l/C by several orders of magnitude. 
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FIG. 8: Spin relaxation rate due to Mt (dashed curve) and M d contribution (dotted curve) to spin-flip transition (see Eq. ( 57 1) 
The sum of these two contributions is plotted by solid curve (R « 1.6 nm, g = 2, T = 0.1 K, V g = hv/150R w 2.3 meV). 



1. Destructive interference 

The remaining and most intriguing interference effect is the one between and Af+ (or between AC and M~). 
These terms are a generated by SO coupling the two lowest states to the excited discrete and the continuous spectrum, 
respectively. From Fig. [6^ and Fig. [5] we see that these contributions at some magnetic field interfere destructively 
leading to a strong increase of the spin relaxation time up to 4 orders of magnitude. Strikingly, such destructive 
interference is robust against a change of parameters, although being most evident when the terms M d and have 
comparable contributions to the spin-flip transitions (compare Figs. [6^1 and[6]3). 

Let us give a physical explanation for this phenomenon. First of all we note that the diagonal elements of the 
electron-phonon coupling (cx g\) (see Eq. (51 1 for details) give the main contribution to the spin- flip transitions with 
respect to the non-diagonal ones (oc 52), since gi ^> gi- As a result, the destructive interference occurs due to diagonal 
elements of I4i_ p h and the elements cx 52 just modulate the strength of the effect, i.e., the depth of the dips in the 
spin relaxation curve. Therefore, in this subsection, we consider diagonal electron-phonon coupling (oc g\) only. 



The terms of the spin- flip matrix clement due to coupling to the first exited subband with m — 1 (see Eq. ( 57 1 ) 
can be written as follows: 



M+ = (AJ)*<* +ifc0jl/a |K 1 -phl^ lfc a 1 i/2> « / Ke^l+ko (0<V, fc0 (0, (58) 

J —OO 1 

1 poo poo poo 

M c = - I ^(A^<^M/2l^i-ph|^ o+ , fco , 1/2 > cx / dCe^< dk{\+y*l + fe (C)^ 0+ , feo (C) (59) 

J kj J —00 J kT 1 ' 
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Note that 3>^ + (C)$ K + k (C) i s a symmetric function of ( with respect to the center of the NT QD (£ = L/2). In 

the dot area (0 < ( < L), it might be approximated by a function cx cos(k (L/2 — £)) with exponential tails in 
the classically forbidden areas (C < and Q > L). In addition to the selection rule Eq. ( 39 1 , it is easy to find that 

( $ K±„fc ,(C)l $ K ±, fc „(C)> = 1 + sgn(m' - l/3)sgn(m - l/3)(-l) n '+™ because $J_± fc ,(C)$ K ±, fcn (0 is either odd or 
even with respect to inversion at £ = L/2. Therefore, Mj~ = at qo = 0, since <I>^ + fc {(/)<& ^+ /- (C) is an asymmetric 
function with respect to (, = L/2 at which it has a node. Thus, & + (0^+ & (0 is found to be well approximated 

by a function sin(fe'(L/2 — £)) defined at < C < L. Now we consider Eq. (59 1. After integration over k, we could 
assume that the dependence of $ >t + k on £ is integrated out, therefore, is a symmetric function of Q with respect 
to ( = L/2, which we approximate by ~icos(k"(L/2 — ()) defined at < ( < L. Using these assumptions, we get the 
following estimations: 

M+ cx f dCe l ^sm(k'(L/2~C)) = fi(q L)+tf 2 (qoL), (60) 
Jo 

M+ cx f dCe l ^icos{k"{L/2-0) = f 3 (q Q L) + if 4 (q L), (61) 
Jo 

, , f L , , „ -fc'cOS^sin^ +qQS - m k^L cos qsL r 

h(q L) = / dCcosg Csin(fc'(L/2-C)) = 2 2 2 H \ 2 -sin^, (62) 

Jo \ K ) ~ % z 



i L , , „ k' cos K± sin ^ - <j sin K± cos ^ o L 
/ 2 (<Zo£) = / dCsin 9o Csin(A : '(L/2-C)) = 2 ? ^ ? cos ^"' ( 63 ) 



, ri , „, , ^ fe" sin K± cos 2^ _ go C os K± sin 2^ g L , . 

= / dCsin 9o Ccos(fc"(i/2-C)) = 2 2 (fc 2 „ )2 J° ? sin ^T' ( 64 ) 

, L , „, , . . -fc" sin ^ cos +g cos^sin^ g £ 

/,«;„£) = -/ dCcos go Ccos(fc"(L/2-C)) = 2 2 ? _™ ? cos ^"- ( 65 ) 
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FIG. 9: (a) Dependence of the real and imaginary parts of matrix elements Mj" and (see Eq. &57l for details) on the ratio 
of the NT length and phonon wavelength qoL (R w 1 .6 n m , I — 100 nm, g — 2). (b) Approximation of the previous dependence 
by fi(qoL) (i = 1, ...,4) (k'L = k" L = 0.7) (see Eqs. p0l-([65|). 



We have plotted the functions fi(qoL) (see Fig.j^ja) in comparison to the real and imaginary parts of M* d (Fig.j9^). 
There is a good agreement between the corresponding functi ons (plott ed with the same line style) except for the 
region of qoL < 2n for f${qoL) and f^(qoL). From Fig. KM Eqs. (62) and (63l, one can see that Re(M^~) and Im(Mj~) 
have zeroes at k! tanq^L/2 = q^tank' L/2 (which are close to q$L = (2n + l)ir, where n = 1, 2, 3 , . . .), in a ddition, 
Re(M^) is zero at q^L = 2nn and Im(M^) is zero at q$L — (2n — From Fig. [9J Eqs. (64) and (65 1, we get 
that Re(M+) and Im(M+) have zeroes at k" tank" L/2 = q ta.nq Q L/2 (which are close to q^L = 2nir), m addition, 
Re(M+) is zero at q L — 2nn and Im(Af+) is zero at q L = (2n — Zeroes of the above functions determine the 
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regions of qoL in which the sign of those functions is constant, namely, 



Re(M+), h(qoL) < 
Re(M+), f 3 (q L) > 
Im(M+), / 2 (q L) > 
Im(M+), / 4 (q L) < 



for 2rwr < qo£ < (2n + l)7r, 
almost for any go-^i 
almost for any q$L > 7T, 
for 2mr < q Q L < (2n + 1)tt. 



(66) 
(67) 
(68) 
(69) 



From these equations we find that for 2nir < go-k ^ (2n+ l)7r the functions Re(M d h ) and Re(M +), as well as, Im(Mj~) 
and Im(Af +) have opposite signs. In other words, due to odd or even symmetry of vector states with respect to the center 
of a NT QD (C, = L/2), the terms of spin-flip transitions Mf and combine in antiphase at 2mr < q$L < (2n+l)7r, 
resulting in destructive interference of those contributions to the spin relaxation rate. Note that ' terms have similar 



behaviour and the same statements hold true for those. Here qo = y V2\hujz — 2r3Ac Urv |/ficsi? is the phonon wave 
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FIG. 10: The dependence of qoL on magnetic field B (R w 1.6 nm, g = 2, Ac Urv = 0.1085 meV). Non-shaded areas stand for 
regions of magnetic field where destructive interference occurs. 



vector of the resonant spin- flip transition. The magnetic-field dependence of go is shown in Fig. 10 (red curve). The 
regions of q^L and B, where the destructive interference is expected, are shown by non-shaded areas. From Figs. [6^, 
[4] a nd [5] we see that the destructive interference dips in the spin relaxation rate are in the defined regions shown in 

Fig. [Tor 

The effect is stronger for smaller V g , when the number of discrete levels is lower (see Fig. |6k). In this case, the 
spin relaxation is predominantly due to coupling to the continuous spectrum. With increasing the voltage V g applied 
to top gates, the number of discrete levels and the spacing between the ground state and the lower bound of the 
continuous spectrum increases. (For instance, in the case of V g = 2.3 meV, there are only two discrete levels, while, 
for V g = 113 meV, there are about 15 quantized levels in each subband.) This decreases (increases) the contribution 
of the continuous (discrete) spectrum to the spin relaxation rate and increases the total spin relaxation rate (compare 
Fig. [5] where V g ss 8.5meV, with Fig. [6^,, where V g « 2.3meV, for the same type of NT QD). Such rich and unexpected 
behavior of the spin relaxation in NT QDs is remarkable and opens up broad perspectives for spintronics in carbon 
nanostructures . 



VII. CONCLUSIONS 



In conclusion, contrary to the common believe that spin-orbit interaction is weak and insignificant in carbon 
materials, we have shown that the situation is actually much richer and that spin-orbit interaction can be very 
important in nanotubes. We have studied spin relaxation and decoherence caused by electron-lattice and spin-orbit 
interaction and predict striking non-monotonic effects induced by magnetic fields B. For particular values of B, 
destructive interference occurs resulting in ultralong spin relaxation times T\ exceeding tens of seconds. For small 
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phonon frequencies to, we find a 1/a/S spin-phonon noise spectrum - a novel dissipation channel for spins in quantum 
dots - which can reduce T± by many orders of magnitude. We show that nanotubes exhibit zero-field level splitting 
caused by spin-orbit interaction. This enables an all-electrical and phase-coherent control of spin - the hallmark of 
spintronics. 

We thank G. Burkard, S. Hani, L. Kouwenhoven, and L. Vandersypen for useful discussions and E. Klinovaya for 
pointing out several typos. We acknowledge support from the Swiss NSF, NCCR Nanoscience, ONR, and JST ICORP. 

APPENDIX A: SPIN SPLITTING AND ELECTRIC-DIPOLE SPIN RESONANCE AT ZERO MAGNETIC 

FIELDS 



From Eq. ( |36[ ), one can find that there is a zero-field splitting between spin-up and spin-down states: 

(£0,0,1/2 - £0,0,-1/2) I s=0 = hvy (-T3V/3R + ALrv/M 2 + k l - hv \J (-T3V/3R - ALrv/M 2 + k o 

« -2r 3 ^ALv, (Al) 
where v = 0, ±1, and where we have taken into account that ko <C 1/R and neglected intervalley mixing. Due 



to the first term in Eq. (34 1, there is spin mixing and, therefore, coupling between corresponding states (see also 



Eq. ( 37 1 ) which allows electric dipole transitions between them. Consider an oscillating electric field (see Fig. Ill: 
E(i) = Ee±smujt, e± is the unit vector perpendicular to the NT axis. An interaction between the electric field and 
an electron in a NT, which leads to electric-dipole transitions, is given by the following operator: 

\e\E „ -i\e\hE d , A „, 

He = cos utF± = — — cosorfsm w— , (A2) 

niQLU moRuj dtp 

where mo is the bare electron mass and P± = —ihsiinpd v /R is the electron momentum along ej_. Here we assume that 
the influence of the lattice potential can be neglected for the estimation of the electric-dipole transitions. Therefore, 



using Eq. p7|, the matrix element of the electric-dipole transitions can be expressed as 

dtp 



, ■ —i\e\hE -r-^ 
(%,0,+1/2 -nB ^0,0,-1/2) = 7; 5— coswt > 



(KY (*„ Zl ,k^i/2\^~' lV Qz\y^M,-i/J 



^(*K b +,* ,+i/al ie i *'^l**+ 1J *,i/a) 



= fhuji cos ut, (A3) 



UR = 2^jS[( A *)* , ( $ *Z 1 ,fc,-l/2l*x b -,fc ,-l/2) + ^( $ ^+,fco,+l/al $ *+ 1 ,fc,l/2>]- 

Here the sum includes summation over the discrete fc„ and integration over the continuous |fe| > 
(\V g \/hv)y/l + 2hv\n rn \/V g . Numerical evaluation leads to the following estimates for the resonance frequency 

u = 2ALrv/fi ~ 33 x 10 10 s" 1 (aLv ~ 0.11 meV) and Rabi frequency u R w 1.6 x 10 5 s" 1 at E = 10 V/ cm 
and V g = 2.3 meV. 

APPENDIX B: VALLEY-ORBIT AND SPIN-ORBIT INTERACTIONS 

In this section, we consider a particle in a NT described by the Hamiltonian ^ with the longitudinal confinement 
([9]) in a parallel magnetic field. The discrete spectrum of a such system is given by 

E m ,n,S ( = i^V Km + fc 2 + S ( huj z , (Bl) 
where K m = (m — + <pab)/R and </?ab = ®ab/®o- Each level is four-fold degenerate (at B = 0) due to valley 



and spin degeneracy. Now, we take SOI into account. For definiteness, we consider only the second term in Eq. (34 1 
which leads to zero-field splitting: 

^ rv - ALvWSc. (B2) 

Moreover, within a minimal model we take intervalley mixing due to non-magnetic impurities or structure defects 
into account. In this case, the intervalley mixing can be described by the following term: 

flk-K' = Ak-k'Ti, (B3) 
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FIG. 11: Scheme of a NT with two electric side gates. An applied ac voltage between the gates creates an oscillating electric 
field perpendicular to the NT axis. Such a setup enables electrically-induced coherent spin manipulation, due to zero-field spin 
splitting in semiconducting NTs. 



where t\ is the Pauli matrix operating on valley-index space. The eigenvalues of the operator Ho + V(£) + H, 
-Z/k-K' (for m = subband) are given by 



E ^ S( = ±hv[v 2 /9R 2 + k 2 n + k 2 l _- K , + (^ AB /R + 2S c k so )' 



2^J^ AB /R+2S c k so )^/9R 2 + k 2 K _ K ,(^/9R 2 + k* n ) 



1/2 



S c tiuiz, 



curv 
SO 



(B4) 



where fcK-K' = Ak-k'/^i k 



so 



Acurv/ftt 1 . The energy spectrum of the lowest electron energy levels and highest 



hole levels described by Eq. (B4| is shown in Fig. 12 (n = 0). In Eq. (B4|, the plus (minus) sign corresponds to 



electron (hole) energy levels. (3=1 for the upper branch of the energy spectrum (blue dashed and red solid curves) 
and (3 = —1 for the lower branch (blue solid and red dashed curves). Using |/t m | ^> k n , we rewrite Eq. (B4| in the 
following way: 



kso) 2 + S c Huj z + O (3vk n R(fi A B) 



(B5) 



From this equation we find that the zero-field splitting of levels is given by fiv^k^ + k 2 ^ ^, and anticrossings (with 

the magnitude 2|A K -K'|) occur at <^ab = — 2Sf&so-R = 0. 

Due to intervalley coupling, electron states of different nonequivalent K points are mixed: 



K 



(0) 
K 



A 



K-K' 



urv J 



(0) 
K' • 



(B6) 



The mixing is maximal at the anticrossing points (at that point, the electron state is just a superposition of those 
corresponding to different K points) and suppressed away from them. Therefore, there is a way to control the 
intervalley mixing for a NT by a magnetic field and it make s NTs attractive for a valley-qubit realization (qubits 
whose quantum states are defined by the valley index) E ? l 38 l 39 l 
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FIG. 12: Magnetic field dependence of energy levels for electrons (blue curves) and holes (red curves). Dashed curves are the 
excited states in a NT QD (R w 1.6 nm, L = 100 nm, g — 2, V a ~ ±8.5 meV (the upper/lower sign is for electrons/holes), 
ALv =0.11 meV, A k _k' = 0.05 meV). 



At anticrossing points, mixing of K and K' valleys is strong and intervalley scattering could occur. Note that 
such scattering is energetically forbidden for first-order processes (with one-phonon scattering). Indeed, scattering 
from K to K' point requires a large change in the electron wave vector (|K — K'| = |K|), while the energy difference 
between the scattering states is small (2Ak-K' < 0.5meV). Phonons in a NT at the K point of the phonon dispersion 
have much higher energy {u>k > 600 cm -1 which correspond to 75 meV)02l and, therefore, first-order single-phonon 
intervalley scattering is forbidden. However, Raman spectroscopy has shown that such an intervalley scattering is 
allowed for photo-excited electrons. Such transitions are attributed to second-order Raman processes by two phonon 
emission or emission of one phonon and elastic scattering on lattice defects (so called D- and G '-bands in the Raman 
spectra) P^l We assume that similar processes could occur in our case, due to spontaneous phonon emission and 
absorption with Hlo\ — hu>i = 2Ak-k' and qi — q2 = K or to emission of a single phonon with q = K and elastic 
scattering on lattice defects, but are nevertheless less probable. 
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